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ABSTRACT 

We explore the assembly history of the M31 bulge within a projected major- 
axis radius of 180" (~680 pc) by studying its stellar populations in Hubble Space 
Telescope WFC3 and ACS observations. Colors formed by comparing near¬ 
ultraviolet vs. optical bands are found to become bluer with increasing major-axis 
radius, which is opposite to that predicted if the sole sources of near-ultraviolet 
light were old extreme horizontal branch stars with a negative radial gradient in 
metallicity. Spectral energy distribution hts require a metal-rich intermediate- 
age stellar population (300 Myr to 1 Gyr old, ^Zq) in addition to the dominant 
old population. The radial gradients in age and metallicity of the old stellar 
population are consistent with those in previous works. For the intermediate-age 
population, we hnd an increase in age with radius and a mass fraction that in¬ 
creases up to 2% at 680 pc away from the center. We exclude contamination from 
the M31 disk and/or halo as the main origin for this population. Our results thus 
suggest that intermediate-age stars exist beyond the central 5" (19 pc) of M31 
and contribute ~1% of the total stellar mass in the bulge. These stars could be 
related to the secular growth of the M31 bulge. 
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Introduction 


As the nearest massive spiral galaxy (~780 kpc, iMcConnachie et al.ll2005l) . the An¬ 
dromeda galaxy (M31) provides an excellent laboratory with which to stndy the constrnction 
of galactic bulges. Compared to the bulge of the Milky Way (MW), the advantages of the 
M31 bulge he in that i) stars are at a single, well-determined distance and ii) it is free of the 
strong, spatially non-uniform line-of-sight extinction against the Galactic bulge. Indeed, the 
M31 bulge is transparent across almost the entire electromag netic spectrum, allowing fo r a 
panchromatic view of its stellar (and interstellar) components (ILi. Wang fc Wakkeiil2009l. for 
a recent review). These advantages largely reduce the uncertainty in deriving the intrinsic 
spectral energy distribution (SED) of the bulge, which is critical to determining the ages 
and metallicities of the underlying stellar populations. 

The stellar populations in the inner bulge of M31 are quite different from those of 


our Galaxy. Around the Galactic cente r, there are a num 


3er o f giant molecular clouds 


few signs of massive stars in the inner bulge of M31. 
near-infrared (NIR) color-magnitude diagrams (GMDs' 

Olsen et ah ( 

of several bub 

2006) investigate the 
je regions, >2' (~450 

pc) awav from M31*, the central super-massive black hole ( 

dress 

er & Richstone 1988: 

Kormendv 1988: Grane et ah 1992: Garcia et ah 2010; 

Li et ah 

2011 

, and hnd that stars 

there are old (>6 Gyr) and metal-rich (>Zr7^) (see also S 

teohens et a 

. 2003: Saraiedini &: Jablonka 


2005h . ISaglia et al.l (120lol) obtain long-slit optical spectra (4787 A—5857 A) in the central 


5' of the M31 bulge, and in most of these spectra hud that multiple spectroscopic line indexes 
can be htted with a single metal-rich stellar population of age around 12 Gyr old. Based 
on the ultraviolet (UV) GMD from Hubble Space Telescope (HST) Faint Object Gamera 
(FOG) observ ations of the c entra l 14"xl4" (53x53 pc) held of M31 in two bands (F175W 
and F275W), iBrown et ahl (119981) suggest that detected UV point sources are not young 
m assive stars, but mainly old evolved low-mass stars, a result conhrmed over a larger area 
by Rosenheld et ah (2012). 


These evolved low-mass stars, i.e., core helium burning extreme horizontal branch 
(EHB) stars and their descendants, have high effective surface temperatures Tg/f > 10^ 
K and can contribute substantially to the UV emission of the M31 bulge. When evolved 
low-mass stars leave the EHB phase, they do not have signihcant stellar envelopes and thus do 
not evolve into the asymptotic giant branch (AGB) phase, and instead either leave the AGB 
phase early (PE-AGB) or fail to enter the phase entirely ( AGB-manque). T hese hot-post HB 
stars (HP-HB stars) are relatively long-lived (~10 Mvr. IO’Gonnelllll999h . In contrast, the 
stars with enough envelopes will complete the AGB phase and then the post-AGB (PAGB) 
phase. PAGB stars are more luminous than the PE-AGB and AGB-manque stars, but be- 
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cause of their short lifetimes (~10^ yr, see review bv lO’Connelllll999l) . PAGE stars contribute 
little to the integrated emission in the UV bands. Indeed, b a sed o n the latest HST Wide 
Fi eld Camera 3 (WFC 3) UVIS observations, iRosenheld et ahl (120121) support the argument 
of Brown et ah ( 1998 1 that most of the detected UV sources in the F275W and F336W 
bands should be HP-HB stars, which can a ccount for only 2% of th e integrated bulge light in 
the F275W band. iBrown et ahl (119981) and IRosenheld et al.l (l2012ll suggest that the remain¬ 
ing extended UV emission arises from the combination of numerous unres olved EHB stars. 


which must be present to explain the HP-HB populations that is observed. IRosenheld et ah 


(120121) also hnd that the fraction of stars passing through this channel increases towards the 
center of the galaxy, though it is never more than 3%. They attribut e this trend, as we ll 
as the positive gradient in the GALEX FUV—NUV color distribution (iThilker et al.ll2005l) . 
wh ich are seen o yer 3 .5 kpc, to the negative radial gradient in the stellar metallicity found 
bv ISaglia et al.l (l2010l) . on the theory that metallicity ahects the strength of stellar winds 
and hence the mass of stellar envelopes and consequently the EHB and HP-HB populations . 
The FUV— NUy color measures the slope of the so-called UV upturn (IThilker et al.l 12005 


0’Connelllll999L and references therein). 


T 

where 


le source o f UV emission becomes more complicated at the very center of M31, 


Lauer et al.l ( 19981) suggest that the FUV peak at the dyna mical center of M31 f ound 


bv iKing et al.l (119951) is slightly extended and is not due to an AGN. iBender et al.l (120051) fur¬ 
ther hnd that the optical spectrum extracted from the central 0.3" held (1.1 pc) is consistent 
with a population of A-type stars or a 200-Myr-old starburst. With observati ons from the 


HST Advanced Gamera for Surveys (AGS) High-Resolution Ghannel (HRG), iLauer et ah 


(I 2 OI 2 I) resolve this central star cluster into twenty near-UV (NUV) bright stars and ar¬ 
gue t hat they can n ot be old stars in the PAGE phase, whose evolutionary duration is 
short. ISaglia et al.l (120101) also show that an old stellar population is not consistent with 
their derived spectroscopic line indexes within the central 5" (19 pc) region and that an 
intermediate-age stellar population is needed. The origin of these intermediate-age stars 
(several hundreds of Myr old) remains an enigma. ISaglia et al.l (120101) suggest that a gas-rich 
minor merger may have triggered a nuclear starburst. If this scenario is valid, it is natural 
to ask whether or not the same event could have triggered star formation in a larger portion 
of the M31 bulge; how much stellar mass this event could have contributed to the bulge and 
what is the event’s role on the evolution of the whole bulge. 


In this paper, we concentrate on studying the stellar populations in the inner bulge of 
M31 with the latest multi-band HST dataset and using our results to constrain the assembly 
history of the M31 bulge. We utilize HST WFG3 and AGS observations in ten hlters from the 
near-UV to near-IR (listed in Table [T]), which cover the central 250" x 250" (~950x950 pc) 
region of M31. Previous photometric observations obtained by ground-based telescopes in the 
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optical/near-infrared bands suffer from confusion and differential extinction due to limited 
and varied resolutions and sensitivities, which hamper the study of the spatial variation 
of underlying stellar populations. Images in our HST dataset have superb resolution and 
sensitivity and cover a broad wavelength range. The inclusion of the near-UV hlters in 
particular is crucial to break the age-metallicity degeneracy in the SED htting. Uncertainties 
in derived ages and metallicities are smaller by more than 60% , when the UV bands are 
inclu ded, than when only optical bands are used, based on studies of 42 globular clusters 
from iKavirai et ahl (120071) . These authors hnd that SED htting using UV photometry also 
outperforms the spectroscopic index method. Further the intermediate-age stellar population 
are expected to have SEDs that peak at 4000 A, making the UV as important constraint. 
Using hts of this full spectrum coverage, we hnd evidence for an intermediate-age stellar 
population. 


The reminder of the paper is organized as follows. We describe the observations and 
data reduction in ^ and the color gradients along the minor axis of the M31 bulge in §3] 
Our method and analysis procedure are presented in ^ and in ^ respectively. We discuss 
our results in §6] and summarize our conclusions in )|71 


2. Observation and Data Reduction 


2.1. HST Multi-Wavelength Dataset 


The maj ority of our dataset comes from the Panch romatic Hubble Andromeda Trea¬ 
sury (THAT; iDalcanton et al.ll2012l: IWilliams et al.l 120141) . As an HST legacy program, this 
survey mapped roughly one third of the M31 disk with the HST WFC3 and ACS in six 
bands, F275W, F336W, F475W, F814W, FllOW and F160W. Except for the FllOW hlter, 
various dither patterns (see Table [T]) were used to remove bad pixels and cosmic-rays, as well 
as to alleviate undersampling. The THAT data include the observations in both the shortest 
and longest wavelength bands (F275W: 2700 A and F160W: 1.5 pm) used in this work, 
giving important leverage needed to constrain extinction and age. F336W and F475W on 
the blue and red sides of the 4000 A break provide a sensitive age measure of the underlying 
stellar populations. The F814W, FllOW and F160W bands, being at longer wavelengths, 
are insensitive to extinction and age and provide a useful color indicator of metallicity. 


We include observations taken in four additional UV/optical bands to increase our 
wavelength coverage. Three of them come from Program GO-12174 (Li et al, in preparation): 
two medium bands (F390M, F547M) and one narrow band (F665N). Their exposure times 
are longer than those of the PHAT program. As a result, signal-to-noise ratios (S/Ns) of these 
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images are comparable to those of the wide hlters used for PHAT. These additional single¬ 
pointing observations covered a limited area of sky. We use the F547M band as an analog 
of the traditional Johnson V band, because of their similar central wavelengths. The fourth 
hlter, F435W, comes fro m observations taken as part of three different programs: Program 
10006, 10760 and 11833 (jWilliams et al.ll2005h . These observations were originally taken to 
monitor the optical counterparts of X-ray novae in the M31 bulge. With various pointing 
directions, rotation angles and exposure times, the F435W observations covered roughly the 
held of interest, except for its southwest corner. The observations were all taken with dither 
patterns. Fig. [1] gives the lo w-resolutio n UV and optical spectrum of the central 10"x20" 
region of M31 obtained from lO’ Connell! fjl999l) . with overlaid the transmission curves of six 
UV and optical bands in our dataset, as well as those of two GALEX bands; FUV and NUV, 
to demonstrate different spectral features covered by these hlters. 


2.2. Data Reduction 

The raw data and calibration hies were downloaded from the Multimission Archive at 
STScI (MAST). The basic calibration of each dithered exposure was made with the latest 
HST pipeline, OPUS version 2010_4, CALACS version 5.1.1 and CALWFC3 version 2.1. The steps 
included identifying bad pixels, bias correction, dark subtraction and hat holding. 

The ‘Astrodrizzle’ task, based on PyRAF0, was used to register individual dithered 
exposures, correct for distortion, mask out defects (including cosmic rays), and combine 
the dithered exposures into pointing images. While ‘Astrodrizzle’ only corrected for the 
distortion of the FllOW images, cosmic-rays in individual FllOW images were identihed 
in the CALWFC3 step, through the use of the sequential readouts ‘MULTIACCUM’ mode. For 
the F275W and F336W observations, in which the WFC3/UVIS chip gap was covered by 
only one dithered exposure, cosmic rays could not be readily identihed. For simplicity, we 
manually removed the regions covered by the gap before mosaicking the images. 



etry of any of the near-UV or optical image relative to the F475W band, using commonly 
detected stars. We aligned the near-IR FllOW or F160W images to the astrometry-corrected 


‘PyRAF’ and ‘Astrodrizzle’ are the product of the Space Telescope Science Institute, which is operated 


by AURA for NASA. 
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F814W images (which were themselves aligned to the F475W images and have more stars 
in common with the near-IR bands). The bias offsets among different pointings in these 
nine hlters were also removed. Finally, we merged the images to form the mosaics in the ten 
bands. These mosaics cover a central 250"x250" (~950x950 pc) field of M31. The ‘PHOT- 
FLAM’ values listed in Table [Dare used to convert the image units from counts (electron/s) 
to flux density (erg cm“^ s“^ A“^). 


Well-characterized uncertainties are crucial for the SED htting. The photometric uncer¬ 
tainty at each pixel in a mosaiced image (a^; the subscript stands for the n*^ hlter) consists 
of three main factors: 1) the Poisson statistical noise of the electronic signal; 2) the error 
in the bias and flat-field calibration; and 3) the systematic uncertainty in the photometry, 
as described by the hts header ‘PHOTFLAM’ keyword. We emp irically characterize the 
factors 1) and 2), using the method described in iDong et ahl (120141) . Briefly, for each pixel, 
we calculated the median and 68% percentile of fluxes in its adjacent 0.65" x 0.65" held as 
the estimates of the local background and its uncertainty. The median ratio of the statisti¬ 
cal uncertainties to the huxes of individual pixels in each of the ten-band mosaics is listed 
in Table [H The WFC3 and ACS systematic calibration error^ are ~2%. We found that 
with the same incoming spectrum, the huxes derived from the SYNPHOT in Chips 1 and 2 
of WFC3/UVIS in the F275W band could diher by up to Therefore, in l|5l we adopt 
the mean value of the F275W huxes from the SYNPHOT for Chips 1 and 2 of WFC3/UVIS, 
including an extra 1.5% systematic error to the F275W band during the SED htting. The 
systematic errors due to the factor 3) are also listed in Table [H In §51 when we present the 
SED htting for a selected region, where the integrated hux uncertainty in each band is the 
square root of the quadratic sum of the empirical statistic errors in all the included pixels in 
this region, together with the systematic errors. 


We corrected for foreground Galactic extinction (E(B—V)=0.062, i.e., Av=0.17. ISchlahv fc Finkbeiner 


20111) for all bands. The relative extinctions {A ^/where A„ or Ap^^pM is t he absolute 
extinc tion in the or F547M band) derived from the MW-type extinction law in iFitzpatrick 
(119991) for an old metal-rich stellar population (12 Gyr ari d solar metallicity) are listed in 
Table [T] As discussed in the Appendix A of iDong et ahl (l2014j) . the relative extinctions 
[An!A p^ 4 Pm) in the ten HST bands are insensitive to the adopted background spectrum. 


^WFC3: http://www.stsci.edu/hst/wfc3/phot^p_lbn and ACS: Bohlin, R., 2012, acs, rept, IB. Private 
communication with STScI helpdesk 

•^The difference could be due to ‘SYNPHOT’, HST/WFG3 STAN Issue 11 
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2.3. Region Selection 


To break the well-known degeneracy between extinction, age and metallicity, we selected 
only low extinction regions for the present study. Fig. [2] presents a map of NIR-to-NUV 
(F160W/F336W) intensity ratio for the inner bulge. This ratio increases with the foreground 
extinction, as well as with the age and metallicity of unresolved stars. Fig. [2] shows that the 
ratio is generally hi gh in regions with hot dust, as traced by the Spitzer/IRAC ‘dust-only’ 8 
/im intensity map of iLi. Wang fc Wakkerl (120091) . This conclusion breaks down in the central 
10" (~38 pc) region, where there is a high ratio of F160W to F336W, but no corresponding 
excess in 8 pm emission. We discuss this region in ( 15.21 A previously known dusty clump 
located at 30" (~120 pc) southeast of M31* (^as outlined by the blue ellipse in Fig. El) 


in CO observations ( 

Melchior & Combes 

Melchior & Combes 

(2013) suggest that 


this dusty clump is located at the far side of the bulge. 


The mean ratio in the South East quadrant of Fig. El is smaller than those in the 
other three quadrants by at least 10%, indicating low foreground extinction. We isolate two 
sectors in this region to avoid the dusty clumps apparent in the flux ratio map and in the 
Spitzer/IRAC 8 pm observations. We use the region between PA = 110° to 180° (green 
sector) for the central 100", while beyond this radius, we use the region between PA=150° 
to 180° (blue sector), to avoid a dusty clump at PA~140° and 160" in projected radius away 
from M31*. According to Li et ah (in preparation), the two-dimensional light distribution 
in the M31 bulge in the ten HST bands can be characterized by elliptical isophotes with a 
PA=50.9° (East of North) and an axis ratio of 0.8. We divide these two chosen sectors into 
elliptical annuli steps of 5" along the major-axis for the SED htting in ( 15.21 We will refer to 
“radius” in this paper as being major axis length of each ellipse. 


3. Color Gradients 

We hrst analyze the radial gradients in UV-optical colors, which give us hints about the 
spatial variation of the properties of underlying stellar populations. Fig. [3] presents the HST 
WFC3/ACS colors obtained from the light intensities in four bands and along the minor-axis 
of the M31 bulge out to 144", corresponding to a major-axis distance of 180" (680 pc) and 
restricted to the two sectors we are using for analysis. In the same plot, we overlay the 
radial surface brightness distribution of the Spitzer/IRAC ‘dust-only’ 8 pm emission. The 
colors tend to become bluer when moving away from M31* and are not correlated with the 
dust distribution. Therefore, the color gradients cannot be solely caused by the decrease 
of extinction and most likely represent a change in the intrinsic properties of the stellar 


















populations. 


Our HST UV-optical color gradients differ from those of GALEX FUV—NUV in iThilker et al. 


(120051) ■ Fig. IHshows the radial color distributions of: F275W—F475W and GALEX/FUV—F475W, 
as well as GALEX FUV—NUV. If the integrated flux in the F275W band was dominated by 
the EHB stars, as suggested by Rosenfleld et al. (2012), the F275W—F475W color, just like 
the FUV—F475W color, can then be used to constrain the ratio of EHB to main-sequence 
turnoff stars. We would expect that the FUV—F475W and F275W—F475W colors follow 
the FUV—NUV color, becoming increasingly red with the galactocentric distance, due to 
the decreasing contribution from EHB stars. But in the fact, the F275W—F475W color 
turns bluer when moving away from M31*; although the FUV—F475W color shows a posi¬ 
tive radial gradient within the central ~50" (190 pc), it turns negative beyond that. These 
different color gradients suggest that the same population of stars that dominates the FUV 
flux (presumed by EHB stars in the old stellar population) probably dominate the FUV 
channel, but do not dominate the intensity in F275W and that instead the variation of the 
F275W—F475W color along the radius may be due to another stellar population (see further 
discussion in l lO.l.ip . 


4. The SED-fitting Method 


We compare the observed and theoretical SEDs to constrain the properties of underlying 
stellar populations. To do this, we minimize the deviation between the observed flux in the 
n*^ Alter, J„, and the Sn predicted by a stellar synthesis model of given stellar age and 
metallicity. 


N 






n=l 




( 1 ) 


where /, An/Ap^ 4 pM and M are the fraction of obscured starlight, the relative ex¬ 
tinction, the absolute extinction and the normalization, respectively, an is the uncertainty 
of Ini as described in §2.21 The sum in this expression is over all available Alters (N<10). 


We utilize ‘Starburst99’ flVazquez fc Leithereiil20051) to calculate the theoretical SEDs. 
Padova stellar synthesis models with full AGB treatments (bu t not the EHB stars and their 
descendants) are adopted. We use the default iKroupal (120021) Initial Mass Function (IMF) 
in St arburst99 (d<F oc M~^dM, F= 1.3 for 0.1 Mq<M<0.5 Mq, 2.3 for 0.5 Mq<M<100 


Mq). IConrov fc van DokkumI (120121) suggest that the IMF of the nuclear region of M31 is 


normal, rather than bottom heavy. While ‘Starburst99’ does not include the recipe for EHB 
stars, the presence of such stars should not affect our results, because their contribution in 
our multi-bands is small (see §3] and ( lO.l.ip . We choose instantaneous burst populations of 
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ages from 100 Myr to 15 Gyr with a step size of 0.01 dex (Myr) and five default metallicities 
(Z=0.02 Zq, 0.2 Zq, 0.4 Zq, 1 Zq and 2.5 Zq). We generate spectra for these ages and 
metallicities, which are then convolved with the hlter transmission curves using the ‘SYNPHOT’ 
to obtain the corresponding fluxes in the ten HST bands. We then linearly interpolate the 
fluxes from these grids to obtain the fluxes of stellar populations with various ages (from 100 
Myr to 15 Gyr) and metallicities (from 0.02 Zq to 2.5 Zq). 


Before we perform the minimization, we choose to £x two parameters. The hrst 
parameter is the fraction of the star light that is obscu red by dust in M31, 0</<l, which 
is degenerate with An (see Fig. 1 of iDong et al.ll2014j) . However, An is insensitive to / in 
regions of low attenuation. This applies to the region that we select for the SED-htting in 
1 12.31 Therefore, we simply £x /=1 {i.e., assuming that extinction exclusively arises from 
the foreground of the bulge). Setting / free in the SED-htting would not change the best-fit 
parameters, but would increase the uncertainties of the output ages and metallicities. The 
second parameter is the relative extinction. An/A f^^^m- We adopt the average extinction 
curve (H„,/Hf 547 m) der ived from five dusty clumps within the central 1' (~230 pc) radius of 
M31 ( Dong et ah 2014 ; Tabled]), which is characterized by i?v=2.4. We also test MW-type 
extinction curve in 1|5] and our results remain unchanged, again because of low foreground 
extinction to the regions of interest. 


5. Results 


5.1. The Presence of Intermediate-age Stars 


As an illustration of our method, we first select a region in the inner bulge of M31 (blue 
box in Fig. [2|) to examine whether or not an intermediate-age stellar population may exist 
beyond the central 5" (19 pc) of the M31 bulge. This region (5"x5", 19 pcxl9 pc) is ~11" 
(42 pc in projection) southeast of M31* and appears dim in the Spitzer/IRAG ‘dust-only’ 8 
/im image, indicating little local foreground extinction. The size of the region is large enough 
for us to obtain the photometry with statistical uncertainties <0.3% in all the ten bands 
and is sufficiently small to neglect any potential variation of the SED along the minor-axis 
of the bulge f 1 15. 21) . 


We hrst ht the ob served SED with one instantaneous stellar population as suggested 
by ISaglia et al.l (120lOl) . This ht has four free parameters: Afmim. age, me tallicity and 
normalization (M). We search for the best ht with MPFIT f Markwardt 20091) . The ht is 
rather poor, especially for the UV bands (see Fig. dj), with x^/d.o.f.=24.6/6=4.1. The age 
and metallicity of the best-ht single stellar population are ~4.4 Gyr and 2.2 Zq, respectively. 
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with Afutm rsj 0.08. The age is much smaller than 12 Gyr as derived bv ISaglia et al.l (120101). If 


we hx the metallicity to the value of [Z/H]=0.16 (1.4 Zq) derived from Fig.l2 of ISaglia et ah 


(120101) at the corresponding radius, the reduced becomes even larger and the htting is 
still poor for the UV bands. 

We thus add another instantaneous stellar population to try to improve the ht. The 
model now includes three more parameters: the age, metallicity and mass fraction of the 
second stellar population. The resultant y^/d.o.f. of the model ht is ~ 4.8/3=1.6, which is a 
signihcant improvement over the ht with the single population. According to the F-test, the 
ht is improved with the null-hypothesis probability p=0.13. The htted extinction is formally 
zero, consistent with the ratio map in Fig. [2l The new model hts the UV bands well (see 
Fig.E]). The age of the old stellar popula tion becomes 14.8 G yi@, while the metallicity is ~1.5 
^0, similar to the values obtained by ISaglia et al.l fj201ol h The second stellar population 
has an age of ~660 Myr and a metallicity of 2 Zq. We estimate the uncertainties of the 
htting parameters via a Monte Garlo simulation. We generate 100 fake SEDs based on the 
best-ht model and the hux measurement errors (assumed to be Gaussian; 1 12.2p . Each of the 
fake SED is rehtted. The 68% percentiles around the best ht value of a parameter are used 
as its Icr uncertainty. Because the SED predicted by the stellar synthesis model is 
a non-linear fnnction of the input parameters, these errors are not symmetric 
to the best-fit valne. In Fig[6l we show the correlation between the age and metallicity 
of the second stellar population. The age is tightly constrained (~0.15 dex), whereas the 
metallicity has a large uncertainty ((5[Z/H]~0.4). This intermediate-age stellar population, 
though accounting for only ~0.9% of the total mass, contributes more than 50%, 20% and 
25% of the integrated intensity in the F275W, F336W and F390M bands, respectively. 

We also perform three tests to check the reliability of our result. First, we randomly 
select several regions of similar sizes and radius in the southeast of the M31 bulge and obtain 
similar ages/metallicities for the two stellar populations. Next, we check for EHB contamina¬ 
tion in the F275W band. We repeat the two component ht excluding F275W band and hnd 
that the ages and the metallicities, as well as the light contribution by the intermediate-age 
stellar population, vary by less than 20%. Finally, we run the SED htting without the UV 
bands (F275W, F336W, F390M). For a single instantaneous stellar population, the y^/d.o.f. 
value reduc es to 3.9/3=1.3. Th e resultant age and metallicity of 9.8 Gyr and ^'Lq are similar 
to those of ISaglia et al.l (120101) . When we utilize two stellar populations, the mass fraction 
of the intermediate-age stars decreases by a factor of ~10, which indicates that the optical 
and IR data alone give little indication for the presence of the intermediate-age stellar pop- 


"^This value is older than the Universe’s age, 13.8 Gyr. However, considering its uncertainty, 0.15 dex, it 
does not contradict the scenario that these stars formed in the early Universe. 
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ulation. Therefore, we conclude that our hnding of the intermediate-age stellar population 
is hrm. 


5.2. Spatial Variation of the Stellar Populations 


We now turn to investigate the radial variation of the stellar populations in the inner 
bulge of M31. We apply the same SED htting as in 1 15.11 for the annuli dehned in 1 12.31 with 
two instantaneous starbursts: one traces the old stellar population and the other represents a 
younger population. Within the sector outlined in green (see Fig.H]), data in all ten bands are 
available for the SED £t, for which the stellar ages and metallicities of the two populations 
are left as free parameters. For the large off-center sector (lOO'To 180", marked blue in 
Fig. |2]), which has data in only seven bands (six PHAT bands -|- F435W), we fix the age 
and metallicity of the old stellar population to the best-fit values obtained for the annulus 
of 95"— 100" in the green sector, considering that from 100" to 180", the prop erties of the 


old stellar population are likely to be constant (see Fig.12 in ISaglia et al.l 120101) . is 

also hxed to zero, which has little effect on the £t. The similar Monte Carlo simulations 
as used in 1 15.11 give the uncertainties of these parameters, which are dominated by the 
systematic photometric errors in the ten bands (the statistical uncertainty in each annulus 
is small, <0.6%). The htting results are listed in Tableland the radial prohles of the htted 
parameters are plotted in Fig. [3 If the two populations are distributed axisymmetrically in 
the bulge, the total mass of the old stellar population is 1.4x10^° Mq and the mass of the 
intermediated age population is 1.5x10® Mq within the central 180" (680 pc). 


6. Discussion 

Using a relatively extinction-free region in the inner bulge of M31, we have characterized 
the radial variations in age and metallicity of the stellar populations. Due to the limitations 
of the available filters, we use only two instantaneous starbursts to fit the observed SED. 
We have found strong evidence for the presence of an intermediate-age stellar population, in 
addition to a plethora of old stars in the inner bulge of the galaxy. 

In this section, we hrst examine the possibility that the extra UV light, which has 
been attributed to the intermediate-age stellar population, may instead come from either 
old metal-rich, evolved low-mass stars in the bulge or main-sequence turnoff stars in the 
projected galactic disk and halo fl lb.ip . We then discuss why the intermediate-age stellar 
population has not been revealed in previous studies f 1 16.21) . Lastly, we explore the potential 











12 


origin of the population and the implications for the star formation history of the M31 bulge 
in 1 16.3I 


6.1. Potential UV sources 

6.1.1. Evolved Low-mass Stars 


We hnd that EHB stars cannot dominate the unresolved emission in the three UV 
bands: F275W, F336W and F390M. The emission peak of the coolest FHB stars with 
temperature 104.15 flo’Connelll 11999 ) is around 2100 A, near the short wavelength 
edge of the transmission curve of the F275W hlter. Furthermore, the intensity observed 
in the broad-band F275W image is contaminated by the emission at optical wavelengths, 
where the M31 bulge is particularly bright, due to a well known red leak in the hlter. 
This contribution from optical light means that only a fraction of the observed 
F275W emission is available to potentially be explained by EHB stars. The 
longer-wavelength F336W and F390M bands sit on the Rayleigh-Jeans tail of 
the EHB SED, making EHB stars very inefficient contributors to these bands. 
Our finding is supported by the data. In ^ we find that the gradients of F275W—F475W 
and FUV—NUV are totally different, suggesting that the FUV band is the only one where the 
EHB flux dominates. Furth er evidence is pre sent in the UV-optical spectrum of the central 
10"x20" region of M31 from O’Connell f 19991) (Fig. [T]). The spectral intensity decreases from 
the optical toward the UV and then increases shortward of 2000 A, presumably due to the 
presence of EHB stars. We use a blackbody to fit the spectrum below 2000 A minimizing 
the potential contamination from low-mass main-sequence turnoff stars and find that log 
(Te//) 4.45, consistent with the value expected for EHB stars. From convolving the UV- 
optical spectrum and the blackbody with the HST WFC3/UVIS transmission curves, we 
find that the EHB contributes at most roughly 23%, 2.9% and 1.2% of the total fluxes in 
the F275W, F336W and F390M bands, respectively, even in these central regions where the 
fraction of EHB stars is expected to be highest. These values decrease with increasing radius, 
because of a decreasing fraction of EHB stars due to the presumably decreasing metallicity. 


We suggest that main-sequence turnoff stars in the intermediate-age stellar population 
found in 1 15.II explain the radial gradients of the near-UV and optical colors. In 1 15.21 we find 
that this relatively young population contributes at least 40% of the emission in F275W (the 
rest is from numerous low-mass main-sequence turnoff stars in the old stellar population) 
in the central 5". This contribution reaches ~62% at 180" away from M31*. The metal- 
licity of this population seems to be constant beyond 30". Therefore, the increasing mass 
fraction of the population at larger radius is expected to cause the systematic blueing of 
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the F275W—F475W color with increasing radius. In the central 50", EHB stars probably 
still dominate the FUV and NUV emission, which explains the positive radial gradient in 
FUV—F475W. However, beyond the central 50", the EHB contribution decreases while the 
FUV and NUV intensities are increasingly contributed by the intermediate-age stellar popu¬ 
lation. Fig. [8] shows the radial distribution of the light contribution in the observed GALEX 
FUV band for the intermediate-age stellar population. Although with large variation due to 
the limitation of our SED htting, there is still a potential trend in this plot; the light contri¬ 
bution by this population increases from 20% near M31* to ~60% at 180". Therefore, the 
emission from main-sequence turnoff stars in that population offsets the metallicity effect of 
the old stellar population on FUV—NUV. As a result, FUV—NUV increases much slowly 
with radius, compared to the central 50". Thus, the intermediate-age stellar population 
detected here is in addition to the EHB stars suggested in previous works. 


6.1.2. Contamination by Disk and Halo Stars 


Considering that the mass fraction of the intermediate-age stellar population increases 

with radius, one might suspect a substantial contamination from stars in the di sk and/or halo _ 

of M3 1. We quantify the con t amin ation based on existing disk and halo models. iKormendv fc Bender 
fllOOOl) and ICourteau et ahl (120 111) have decomposed the light contributions from the disk, 
bulge and halo of M31 through htting the surface brightness distributions in the U ad J 
bands, respectively. Their decomposition gives the half-light radius of the bulge and the 
scale length of the expo nential disk as 0.8 k pc (210") and 5 kpc (1320") (see also the re¬ 
cent characterization in iDorman et ahl 1201311 . The scale length of the disk is substantially 
larger than the size of our regions. With the small held-of-view considered here, their abso¬ 
lute contamination should be approximately uniform. Fig. [HI compares the model-predicted 
contribution in the V band from the bulge and the disk and our measured contribution 
from intermediate-age stellar population. We hud that within the central ~100", the disk 



then expect that they contribute < 3% of the flux of the intermediate-age stellar population 


in the V band. Halo stars are also too metal-poor to explain the intermediate-age stars. 
Most importantly, the observed mass surface density distribution of the intermediate-age 
stellar population (Fig. ITOl) is very much different from a flat distribution as may be ex¬ 
pected from the disk and halo contamination. Therefore, we conclude that the bulk of the 
intermediate-age stellar population indeed arises from the M31 bulge. 
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6.2. Comparison with Previous Work 


Previous studies of the stellar populations in the inner bulge of M31 use op tical and/or 


near-I R data which are not particularly sensitive to intermediate-age stars. lOlsen et ah 


(120061) analyze the images taken by Gemini/NIRI and HST/NICMOS observations of the 
M31 bulge. Even with the high angular resolution provided by the the Gemini/NIRI with 
Adaptive Optics (~0.09") and HST/NIGMOS (~0.185"), they detect only individu al stars 
at the tip of the red giant branch (RGB) with ages >1 Gyr. ISaglia et al.l (120101) ht the 
ground-based long-slit optical spectra within the central 5" (19 pc) radius with two stellar 
populations. They hnd an intermediate-age stellar population with age ~600 Myr, which 
contributes less than 10% of the total stellar mass, with the rest from the old stellar popu¬ 
lation (~8 Gyr), which is consistent with our result. At larger radii, they only use one old 
stellar population to infer the presence of a negative radial gradient in metallicity, as well as 
the stellar age of ~12 Gyr, a trend we also see from the broad band hts. They then compare 
their Lick index analysis results with the u-g color of the SDSS data and conclude that they 
are consistent. 

The non-detection of the intermediate-age stellar population in bulge regions beyond 
the central 5" in previous works is likely due to their detection limit. According to the 
Padova stellar evolutionary tracks (IBressan et al.N20121) . for a stellar population with Zq 
and age=0.7— 1 Gyr, stars with masses >3 Mq have died, while stars with 2.1—2.4 Mq 
should have just left the main-sequence. The F275W and F336W magnitudes of these latter 
stars (>25 mag) are still fainter than the 90 % completeness limit of the PHAT survey in 


the UV band (~24.5 at F275W and F336W, iRosenheld et al 


explain why they are not detected in the study by ISaelia et al 


20121). Their age can also 
( 2010 ). In the Lick index 


method, the H/5 4861 A absorption line is a critical age discriminator, because it is strong 
in A-type stars of age around few hundred yea rs. Under-subtracting the H/5 emission line 
from the nuclear gas spiral surrounding M31* (iLi. Wang fc Wakkerl 120091) may have led to 


an ov erestimate of the stellar age based on the Lick index method, although LSagha et al. 


(I 2 OIOI) mention that they estimate the properties of the underlying stellar population and 
emission lines, simultaneously. Because of the strong H/5 absorption line, the existence of 
~300—600 Myr old stars near M31* is well recognized. Away from M31*, the age of the 
intermediate-age stellar population detected here is large enough that the corresponding 
decrease in the equivalent width of the H/5 line makes it hard to detect the population in 
the integrated spectrum, which is dominated by the old stellar population. In ( 15.11 we have 
shown that without the three UV bands, the observed SED could indeed be htted reasonabl y 
well with a single stellar population, reaching the same results as those of ISaglia et al.l (120101) . 
Therefore, including the UV bands into the SED htting is critical to the identihcation of the 
intermediate-age stellar population. As a result, we are more sensitive to this population 
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than any previous study. 


6.3. The Building History of the M31 Bulge 


The exis tence of metal-rich in termediate-age stars in the M31 inner bulge is not a unique 
phenomenon. iBensbv et ahl (l2013l) present a study of the high-resolution spectra of 58 dwarfs 
and supergiants in the Galactic bulge. Metal-poor stars are shown to have relatively old ages 
(10—12 Gyr), whereas metal-rich ones ([Fe/H]>-0.1) have a wide range of ages from 2 to 12 
Gyr and at leas t 5% of them are young er than 5 Gyr. Based on a study of eight nearby 
spiral galaxies, iMacArthur et al.l fl2009f) hnd positive radial gradients of luminosity-weight 
age, which indicate the existence of young stellar populations in the galactic nuclei. They 
suggest that relatively young (< 1 Gyr) stellar populations contribute as much as 70% of the 
optical emission from the galatic bulges, although their mass contribution are small, <20%. 

The radial gradients in the age and metallicity of stellar populations provide insights 
into not only their origins, but the formation history of the M31 bulge as well. The stel¬ 
lar bulges of spiral galaxies may be classihed into three types: 1) ‘classical’, 2) ‘disk-like’ 


and 3) boxy/peanut-shaped flKormendv &: KennicuttI l2004 lAthanassoulal l2005l) . Glassi- 
cal bulges have Sersic index >2 and have properties similar to elliptical galaxies. They 
also fall into the fundamental plane of elliptic al galaxies and are supported primarily by 
random motion flKormendv fc KennicuttN2004l) . Such spheroids are suggested to be con¬ 
structed through monolithic co llapses of primordial gas clouds and/or major mergers of 
galaxies (IMacArthur et al.ll2009l. and reference therein). In the monolithic collapse scenario, 
inner parts of spheroids experience an intense star formation in the early Universe over a short 
time interval. Heavy elements released by this activity are locked in the galactic nucleus, 
because of its deep potential well. As a result, stars formed in subsequent star formation 
tend to be metal-rich. Meanwhile, metal-poor intergalactic gas could keep falling into the 
galaxy, which may trigger new star formation in outskirts and produce younger and metal- 
poorer stars. Therefore, this scenario predicts negative age and metallicity gradients. After 
a major merger, galactic disk could be disturbed, or even destroyed to form galactic bulge. 
The age and metallicity distributions in the new spheroid are smoothed and become flat. 
‘Disk-like’ bulges with Sersic index <2 are offen identihed as ‘pseu dobulges’. Boxy/peanut 


bulge s are generally interpreted as edge-on barred pseudobulges ( Bureau fc Athanassoula 
2 OO 5 II . A pseudobulge is likely built through ‘secular evolution’ flKormendv fc Kennicutt 


20041) . The disk instability introduces a rotating bar, which buckles and heats the disk verti¬ 
cally to increase its scale height. The bar can also transport the disk material into the inner 
region, leading to star formation and thus the growth of the bulge. This kind of bulge still 
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keeps the memory of the galactic disk, so that they are primarily rotationally supported. 
Because it takes a long time for materials to spiral into the inner parts of such a bulge, a 
negative age gradient and a uniform metallicity are expected. 


The M31 bulge is classihed as a classical bulge, because of its Sersic index. 2.2 fICourteau et ah 


201 


and 


), although this val ue is close to the boundary to a pseudobulge flFisher fc DrorvI 120101) 


Beaton et al.l (120071) hnds boxy structures. The large velocity dispersion, ~160 km/s fjSaglia et al 


201ol) . also supports that M31 has a classical bulge. While the monolithic collapse scenario 
coul d explain the age a nd metallicity gradients of the old stellar population found in 1 15.21 (see 
also lSaglia et al.ll2010|) . the supersolar metallicity indicates subsequent major merger(s). But 
the presence of the dynamic fragile galactic disk indicates that no recent major merger has 
happened in M31. Therefore, the high stellar metallicity and its negative gradient in the 
M31 bulge suggest that the early M31 experienced several major mergers, which triggered 
the formation of stars with enhanced metallicity. After that, accreted intergalactic materials 
have reduced the metallicity of stars in the onter parts of the bulge. 

On the other hand, the intermediate-age stars may be due to secular evolution, instead 
of mergers. M31 has likely experienced minor mergers recently, with its satellite galaxies, as 


indic ated in stellar number density maps produced by wide-held imaging surveys fllbata et al. 
20 14 and references therein). A potential head-on collision between M31 and M32 about 


2010; 

Lauer et al. 

2012 

), as well as the 10 kpc star forming ring ( 

Block et al. 


20061) : although recently iDierickx et a 


to produce the ring and iLewis et al 


(120141) suggest that head-on collision is not required 


( 2015[) hnd that the 10 kpc ring is long-l ived and 


stationary, which cannot be due to a purely collisional origin. iBlock et al.l (120061) snggest 
that the M32 just passed the inner bulge of M31 recently, which triggered the star formation 
in its nucleus. This scenario predicts that the age of the intermediate-age stellar population 
should be constant in the inner bulge of M31, which is inconsistent with the positive radial age 
gradient. Therefore, although we cannot exclude the possibility that this fly-by interaction 
produced the young star cluster in the central 5", it cannot explain the intermediate-age stars 
further beyond. Most likely, the intermediate-age stellar population was formed in the secular 


evolut i on scenario. The pres ence of a bar in the M31 bulge is reported by iBeaton et al. 


(120071) . ICourteau et al.l (l201ll) also suggest that the rat io of the bulge-t o-disk scale lengths of 
M31 is ~0.2, a prediction of secular evolution models (1Courteaulll9961) . which also naturally 
explains the age gradient and ronghly constant metallicit y of the intermed iate-age stellar 
population. If we assnme a 5% star formation efficiency (lEvans et al.l 120091) and that disk 
gas was smoothly transported into the bulge during the last 1 Gyr, an infall rate, 3.4 MQ/yr, 
is expected. This value falls within within the range predicted by the num erical simulations 
of the evolution of the galactic disks (0.1-10 MQ/yr, iMinchev et al.l|20121) . 
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The pseudobulge and classical bulge could co-exist IjErwin et al.ll2015h . iFisher &: Drorvi 


(120 lOl) find that there are a group of pseudobulges (defined by their low Sersic index and mor¬ 
phology) with very low specific star formation rates, just like the classical bul ges. Adding a 
pseud obulge into the classical bulge could also efficiently reduce its Sersic index flFisher fc Drorv 
20101 ). This may explain the low Sersic index of the M31 bulge, compared to the other clas¬ 
sical bulges. 

Apparently, the suggested secular evolution contributes little to the overall mass of the 
M31 bulge. Assuming secular evolution is responsible for the entire intermediate-age stellar 
population detected, the mass fraction of intermediate-age st ars derived from our SEP fitting 
is indeed small, ~1%. This is consistent with the result in iMacArthur et ahl (120091) . which 
suggest that bulge growth via secular processes generally contributes little to the stellar mass 
budget in their eight nearby spiral galaxies. Therefore, the M31 bulge is similar to the MW 
bulge, in which the dominant majority of old stars is only contaminated by relatively few 
young stars. 


7. Summary 

We have studied stellar populations in the inner bulge of M31 with multi-wavelength 
observations taken with ffST WFC3/ACS. The broad wavelength range coverage of our data 
set from near-ultraviolet to near-infrared enables us to decompose multiple stellar popula¬ 
tions in regions beyond the central 5". We summarize our results below: 

• In the near-ultraviolet to near-infrared range, the light becomes increasingly blue with 
the galactocentric distance in the bulge. This trend is in sharp contrast to the positive 
radial gradient seen in the GALEX FUV—NUV color. The FUV—F475W radial color 
gradient changes at about 50" major-axis radius. These trends cannot be explained by 
the possible presence of old post-EHB stars with a negative radial metallicity gradient. 

• We have found that in addition to the known old stellar population, an intermediate- 
age stellar population is most likely needed to explain the observed SED and its radial 
change across the inner bulge, especially in the three UV bands. Although the mass 
fraction of this new population is low (0.2—2%), its emission is significant (e.g., > 40% 
in the F275W band). 

• The mass surface density of intermediate-age stars decreases with the radius. This 
trend rules out the possibility that they represent the projected galactic disk contribu¬ 
tion of M31. We find that this contribution accounts for at most 50% of the V band 
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intensity in the central 100" region. Meanwhile, the metal-poor M31 halo stars can 
only contribnte at most extra 3% the V band intensity in the same region. 

• The age (~300 Myr—1 Gyr) and mass fraction of the intermediate-age stellar pop¬ 
ulation increase with the galactocentric distance from the center of M31, while the 
metallicity is ronghly constant (^Zq). 

• Because the intermediate-age stars are older than the young star cluster (age ~200 
Myr) in the central 5", their formation is unlikely due to the possible head-on colli¬ 
sion between M31 and M32. We propose that the population represents the secular 
evolution of the inner bulge via star formation from the inflow gas, probably induced 
by the presence of the bar of the galaxy. This secular evolution could also explain the 
positive radial age gradient of the intermediate-age stellar population. 

• The intermediate-age population contributes only ~1% of the stellar mass in the central 
180", indicating that secular growth plays an insignihcant role in building the M31 
bulge. 


The rad ial age and rnetallic ity prohles of the old stellar population are consistent with 
those of ISaglia et al.l (120101) : its age does not show any signihcant variation, while its 
metallicity changes from ~ 2Zq to ~ Zq in the held considered here. This suggests 
that major mergers taking place in the early Universe produced the majority of old 
metal-rich stars. Subsequent star formation from the accretion of metal-poor materials 
may have resulted in the negative metallicity gradient observed. 
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Table 1. HST Multi-Wavelength Observations 


Filter 

Detector 

Pivot A 
(A) 

Program 

ID 

Number of 
dithered exposures 

Exposure 
time (s) 

PHOTFLAM 
ergs cm~^ s~^ A~^ 

Systematic error 
of PHOTFLAM‘^ 

Median of 
a/I (%)>> 

An/Ap54jM 
Milky Way^^ 

An/Ap54jM 
M31 bulge‘s 

F275W 

WFC3 /UVIS 

2704 

12058 

2 

925 

3.3010e-18 

3.5% 

13.1 

1.94 

2.80 

F336W 

WFC3 /UVIS 

3355 

12058 

2 

1250 

1.3129e-18 

2% 

3.7 

1.64 

1.92 

F390M 

WFC3 /UVIS 

3897 

12174 

3 

2700 

2.5171e-18 

2% 

3.1 

1.47 

1.58 

F435W 

ACS/WFC 

4319 

10006 

4 

2200 






F435W 

ACS/WFC 

4319 

10760 

8 

4360 

3.1840e-19 

2% 

3.6 

1.31 

1.37 

F435W 

ACS/WFC 

4319 

11833 

8 

4360 






F475W 

ACS/WFC 

4747 

12058 

5 

1900 

1.8210e-19 

2% 

4.1 

1.18 

1.21 

F547M 

WFC3 /UVIS 

5447 

12174 

3 

2700 

4.6321e-19 

2% 

4.1 

1.0 

1.0 

F665N 

WFC3 /UVIS 

6656 

12174 

3 

2700 

1.9943e-18 

2% 

5.0 

0.75 

0.73 

F814W 

ACS/WFC 

8057 

12058 

4 

1715 

7.0332e-20 

2% 

7.8 

0.55 

0.61 

FllOW 

WFC3 /IR 

11534 

12058 

1 

699 

1.5274e-20 

2% 

12.7 

0.31 

0.46 

F160W 

WFC3 /IR 

15369 

12058 

4 

1600 

1.9276e-20 

2% 

14.7 

0.19 

0.43 


Note. — a) The ‘PHOTFLAM’ information are from: http://www.stsci.edu/hst/wfc3/phot_zp_lbn (WFC3) and http://www.stsci.edu/hst/acs/analysis/zeropoints/ 7 ^tablestart 
(ACS), b) a and ‘F are the uncertainty and inten sity of each pixel, c) The relative extinction An j for the MW-type dust, d) The average relative extinction An !Ap^^-ji^ 
of five dusty clumps in the M31 bulge, derived bviDong et all ^201411 ■ 
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SED fitting for the southwest region of the M31 bulge 
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i-9o_0 81 

1-9-0.5 

1 -^- 0 .7 

1.97 
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0 4+0-3 

9-"1-0.2 
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1.38 
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U.yy_Q gg 

0 c:+0.4 
9.0 _q 2 

1-1-0.4 

1.52 

110-115 

11.8 
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1.10 

39.815-? 
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11.8 
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36.915:5 
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1 oq + 1-22 
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1.37 
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11.8 

1.10 

36.11?-? 
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0 5+0.5 
9.0_q 2 

1 4+0-0 

1-^-1.1 

1.30 
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11.8 

1.10 
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8I1I55? 

0.8215:5? 

9-9_o.3 

1 q+0-2 

1-0-1.4 

1.27 
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11.8 
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34 7+^'^ 

< _o.9 

7631555 

0.9015:5? 

9-9_o.3 

1 7+1-0 

i-'-i.o 

1.19 

150-155 

11.8 

1.10 

33 9+^-^ 

oo.y_o_3 

81115?I 

0-9215:1 

0 7+0.5 

9- '-0.3 

1 9+1-8 

1 -^- 0 .7 

1.14 

155-160 

11.8 

1.10 

27 6^^'^ 

8I1I55? 

0-8515:55 

0 r+0.4 

9-0_o.3 

1-915:5 

1.03 
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11.8 

1.10 

21 
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Note. — The superscript and subscript represent the 68% confidence range. Beyond the 100^^ there are no F390M, F547M 
and F665N observations. In order to have enough freedom for the least fitting, we freeze the age and metallicity of the old 
stellar population to the values at the annulus 95''-100''. 
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Fig. 1.— The spect rum of the central 10"x20" region of M31 (solid line) obtained 
from lO’Connelll fll999l) . The units of y-axis are Vega magnitude. The UV part of the 
spectrum (below 3200 A) is from the International Ultraviolet Explorer, while the part 
above 3200 A is from a ground-based telescope. The dotted line represents the EHB contri¬ 
bution inferred from a black-body fit to the spectrum below 2000 A. The dashed lines are the 
transmission curves of the detectors: GALEX (blue), ffST/WFC3 (green) and ffST/ACS 
(red) (from left to right, GALEX/FUV, GALEX/NUV, F275W, F336W, F390M, F435W, 
F475W, F547M) 
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3.9 




Right ascension 


Fig. 2.— Map of the F160W to F 336W intensity ratio, wit h overlaid contours from 
Spitzer/IRAC 8 /rm dust-only image (ILi. Wang fc Wakkerl 120091 ). The blue circle outlines 
the central 10" region. The blue box marks the region used in §5.11 to constrain the stellar 
population in the off-center M31 bulge. The green and blue sectors are used to study the stel¬ 
lar population variation along the minor-axis of the M31 bulge (see 1 15.2p . Three foreground 
stars and globular clusters (marked with small green circles) are removed. The blue ellipse 
encompasses one dusty clump at 30" (he., ~120 pc) southeast of M31*. The black spot to 
the southeast is the death star (very low sensitivity) feature in the WFC3 /IR detector. The 
two black strips are excluded regions, which were covered by only one dithered exposure in 
the F275W or F336W bands and do not allow for cosmic-ray removal. 
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Fig. 3.— The radial distributions of four near-UV to near-IR colors along the minor axis of 
the M31 bulge (see 1 15.2p . The color uncertainties are smaller than the size of the symbols. 
The distribution of the average surface brightness of Spitzer/IRAC ‘dust-only’ 8 /rm image 
in the same sky (black curve) is rescaled for comparison. 





color+constant 


27 


Major-axis radius (pc) 

100 200 300 400 500 600 



Major-axis radius (arcsecond) 


Fig. 4.— The radial distributions of the GALEX/FUV—F475W, F275W—F475W and 
GALEX FUV-NUV colors. 
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Bands Bands 

Fig. 5.— The x values from the modeling with single stellar population (left panel) or 
intermediate-age+old stellar populations (right panel). 
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[Z/H] 


Fig. 6 .— The relationship between the metallicity ([Z/H]=logio(Z/Z 0 )) and age of the 
intermediage-age stellar population, as inferred from the two-population SED htting for a 
randomly selected region in the M31 bulge (blue box in Fig. |2]). The contours are at the 
68 , 90 and 95% conhdence levels. The ‘plus’ symbol marks the best-£t metallicity and age. 

Limited by the stellar synthesis model, our contours are not extended to [Z/H] 
> 0.4. 
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Fig. 7.— The radial distributions of the stellar age, metallicity, mass and light fraction 
for the young (intermediate-age) and old stellar population, as marked in individual panels. 
The bottom right panel shows the light fraction for the young (intermediate-age) stellar 
population in the integrated light in the F275W and F814W bands. In each panel, the grey 
area represents the uncertainty range at the 68% conhdence level, while the error bars are 
the uncertainties of the means at individual radial intervals. Due to the limit of data 
points (see Section 15.2h . beyond 100", we fix the age and metallicity of the old 
stellar population to the corresponding values in the annulus of 95" to 100". 
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Fig. 8.— The radial distribution of the light contribution for the i ntermediate-age ste llar 
population (FUVj„t) in the observed GALEX FUV band (FUVga^ex, Thilker et ah 2005 1. 
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Fig. 9.— Comparison between various V-band/F547M intensity fractions. The black solid 
and dotted lines repre sent the total fractions cont ributed by the bulge and disk of M31 in the 
C-band, respectively (iKormendv fc Bender! 119991) . The square and error signs represent the 
fraction from the intermediate-age stellar population and its 68% uncertainty in the observed 
F547M intensity as estimated from our SED fitting. Because there is no F547M data beyond 
100", we show only the curves within this region. The red dashed line is for the fractional 
contribution of the disk to the intermediate-age stellar population (he., the squares divided 
by the black dotted line). 
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Fig. 10.— The mass surface density of the intermediate-age stellar population as a function 
of the major-axis radius. 




